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SUMMARY 

Some details of the physical mechanism of the combustion of mag- 
nesiiom ribbon were investigated, and burning times for magnesium ribbons 
were determined with various mixtures of oxygen in argon, nitrogen, 
helium, and argon - water vapor. Evidence from this investigation indi- 
cates that a gas-phase rather than a surface reaction occurs, which 
bears some analogy to the combustion of liquid-fuel drops. The experi- 
mental burning times are compared with those calculated by using a 
vapor -phase heat- and mass-transfer mechanism. The infinitesimally thin 
flame front concept successfully applied to fuel drops was abandoned in 
favor of a high-tenperature reaction zone of finite thickness. Agree- 
ment between experimental and calculated values for burning time in terms 
of oxygen concentration for the argon- and nitrogen-inerted systems was 
quite satisfactory on a relative basis; actual numerical values agreed 
within a factor of 2 to 3. Agreement was less satisfactory for low- 
density helium-inerted atmospheres, although the changes observed ex- 
perimentally were in the direction indicated by the calculation. 

For the geometry employed, observed burning times (approximately 
proportional to reciprocal average burning rates) decreased by a factor 
of approximately 10 with increasing oxygen concentration over the range 
17 to 100 percent by volume. The times were slightly longer with argon 
mixtures than with nitrogen mixtures and appreciably shorter with helium 
than with argon or nitrogen. Burning times decreased by as much as one 
third with increasing moisture content; the greatest changes occurred 
for small amounts of water vapor and for long burning times. 


INTRODUCTION 

In spite of the expanded efforts on combustion research in recent 
years, the burning of metallic solids has not received much attention. 
While solid carbon has been the subject of extensive investigations, 
which clearly indicate the significance of surface reaction in the utili- 
zation of nonvolatile carbon, there is relatively little information 
available on the detailed mechanism of the combustion of metals. A 
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possible reason for this lack of information is that previous investi- 
gators may have generally assumed that the burning of metal solids re- 
sembles that of liquid hydrocarbons. 

Work on the carbon problem has ranged from the isolation of the 
surface reactions (refs. 1 to 3) and the burning of carbon rods and 
single carbon spheres in air streams (refs. 4 to 6) to complex studies 
of fuel beds (ref. l) , coal-gasification projects (ref. 8), and the 
powdered-coal gas turbine (ref. 9). The utilization of solid explosives, 
primarily as ammunition (ref. lO) , is another aspect of solid combustion 
which has received much attention. A project has also been conducted in 
the attainment of high temperature by burning metals (ref. ll) . 

In a somewhat different vein, dust explosions, first of coal, flour, 
and sugar, and more recently of metallic dusts in powder metallurgy, have 
attracted attention from the point of view of safety; ignition energies, 
explosive limits, pressure rises, and rates of pressure rise have been 
studied (ref. 12). Closer to the subject is work on the propagation of 
flame through dust clouds: the experimental flame speed of aluminum 
powder dust together with the effects of concentration and particle size 
(ref. 13), spectral studies of metal dust flames (ref. 14), and work with 
metal dusts in Bunsen biorners including the determination of flame tem- 
peratures (refs. 15 and 16). 

Of much general interest, but little direct significance, is the 
considerable volume of work on the high-temperature slow oxidation of 
metals; that is, studies of the growth of oxide films at temperatures 
below those at which ignition occurs. These investigations (refs. 17 
and 18), covering a wide range of metals, definitely involve metallic 
surface phenomena. Of more specific interest are studies of the burning 
of single liquid fuel drops (refs. 19 to 2l) , and of liquid fuel from 
burners of various geometries (ref. 22). In such cases, where vapor proc- 
esses are known to occur, successful techniques for the treatment of 
single diffusion flames have been developed. A s-ummary of the influence 
of inert gases on homogeneous gaseous flame phenomena has been prepared 
(ref. 23). (The inert gases used in the experiments described in this 
report are included in ref. 23.) 

The purpose of this investigation was to study the gross physical 
mechanism of the combustion of magnesium metal in order to learn whether 
a vapor -phase or siorface reaction is involved; further, to examine the 
effect of the composition of the oxidizing atmosphere, both as a matter 
of general interest and in the hope that such studies might be informa- 
tive as to the details of the physical and/or chemical processes in- 
volved. In particular, magnesium ribbon was burned in 17- to 100-percent 
oxygen mixtures with argon, nitrogen, helium, and argon - water vapor as 
diluents. Data are presented in terms of burning times for a strip of 
ribbon about 4.6 centimeters long; such data bear an approximately in- 
verse relation to average burning rates. The experimental data are 
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compared with average burning times calculated from hypotheses involving 
diffusion and heat-transfer relations; the effects of moisture are less 
satisfactorily explained. The occurrence of a vapor-phase combustion 
process is strongly indicated on the basis of (l) the experimental ob- 
servations and (2) the comparison between experimental and calculated 
burning times. 


w APPARATUS 

cnd 

Box . - All experiments were carried out in a cubical box with an 
inner dimension of 7 inches^ constructed of l/2-inch transparent plastic 
(fig. 1). A rectangular access port in the top was covered by a metal 
plate with a sheet of heavy gasket material attached; the cover served 
as an adequate seal under vacuinn conditions and as a vent against pres- 
s-ures above atmospheric. A rubber stopper in the bottom supported two 
electrodes and a tube for evacuating and filling the box. The tube was 
connected to the gas source, and through a filter to a pressure gage 
and the laboratory vacuum system. The magnesium ribbon was supported by 
a metal bracket with a gap of about 4.6 centimeters and held in place by 
small metal weights. 

Ignition . - The electrodes supported a heater in the shape of a 
squared "U” about 1.5 centimeters across and formed from about 8 centi- 
meters of number 16 Nichrome wire. To produce ignition, a thread, con- 
sisting of one third of a 7.5-centimeter magnesium ribbon, was wrapped 
tightly about the cross bar of the heater and placed just under the 
ribbon; 4 volts across the heater fired the thread, and the thread in 
tiorn ignited the ribbon. 

Magnesium ribbon . - The magnesiimi ribbon measured 0.015 by 0.31 
centimeter and weighed 0.0074 gram per centimeter (about 10 percent less 
than weight calculated from dimensions). No special precautions were 
taken to remove oxide film or prevent oxidation, but the ribbons were 
kept in closed containers and retained their luster. 

Gases . - Mixtures of 17.6, 21.5, 24.1, 34.2, and 43.7 volume per- 
cent oxygen in argon were obtained commercially and the analyses checked. 
The 18.2, 32.0, and 70.6 percent oxygen in argon and the 18.0 and 53.0 
percent oxygen in helium were mixed in the laboratory from cylinders of 
the individml gases in a 35-liter tank to a total pressure of 5 at- 
mospheres. Compressed air and mixtures of 34.9, 44.5, and 70.0 percent 
oxygen in nitrogen were also obtained commercially. 

Water vapor . - Moisture was added to the 21.5 and 43.7 percent oxy- 
gen in argon mixtures by flowing the gas mixtirre over water in a simple 
train; the amount of water taken up was varied by adjusting the ratio of 
gas flowing through the train and through a dry bypass. 
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The concentration of water in the mixture was determined in terms 
of the dew-point temperature; with adequate conditioning of the instru- 
ment^ readings obtained before and after runs of several hours checked 
to within about 5 percent of the amount of water present. 

Cameras . - Motion pictures were taken at about 16 frames per second 
with spring-driven cameras calibrated for speed and constancy of speed. 
Detailed still pictures were obtained with a view camera. 


PROCEDURE 
Obtaining Data 

The procedure used for obtaining data was as follows: The thread 

and the ribbon were placed in position. The box was twice evacuated to 
a pressure of 25 to 40 millimeters and filled with the selected mixture. 
In addition, for mixtures containing moisture, the box was purged for 8 
to 10 minutes (three to four times the period required to fill the box 
at the flow rates used in the train). The only effects of procedure 
were found for mixtures containing water vapor. External electrical and 
photographic adjustments were made while the box was being filled. With 
all preparations completed, the heater current was turned on and then 
shut off when the thread ignited. The motion-picture camera was started 
at the first glow of the Richrome; the still camera was snapped manually. 


Reducing Data 

In the reduction of data, all quantitative results were derived 
from the motion pictures. Selected and numbered frames are shown in 
figure 2 for a test with .a burning time of about 3 seconds. 

Burning time . - The time elapsed between the start and the finish 
of the burning of the magnesium ribbon spanning the 4.6-centimeter gap 
in the metal bracket is defined as the burning time. No corrections 
were made on the burning-time measurements for the action of the igniter 
or for the quenching effect of the metal bracket (burning approached the 
bracket somewhat more closely for shorter burning times). The start was 
taken as the first frame showing a flame which within several frames was 
observed to be spreading (extending horizontally). For example, in fig- 
ure 2 no spreading occurs between frames -10, -5, and 0, but by frame 5 
spreading is underway; therefore, 0 is taken as the start (in practice, 
with a viewer, this procedure is capable of finer gradations). Frame 30 
shows fully developed combustion. By frame 37 the flame is starting to 
subside. An appreciably decreased luminous area is shown in frame 45, 
and frame 50 indicates the final stages of vigorous burning. Frame 53 
shows a luminosity which thereafter changes relatively slowly with time 
(it persists considerably beyond frame 59). Therefore, frame 53, as an 
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intermediate condition between vigorous burning and residual luminosity, 
is taken as the end. On the whole, the reproducibility of this measur- 
ing procedure is apparently quite adequate (the spread of the data is 
shown in fig. 5); the data are of such a nature that a systematic error 
of a few frames, either as a constant or as a percentage error, will have 
little effect in the final comparisons. 

Apparent flame diameter . - The distance from the bottom of the lumi- 
nous zone to the point above the ribbon where the luminous intensity 
starts to decrease is the apparent flame diameter (frame 30, fig. 2), 
an average over a number of points on each frame and over a number of 
frames. The numbers are at best relative on a rather indefinite scale 
because of constant exposure times and differences in energy-release 
rates producing some overexposure, and because of the difficulties of 
visually estimating intensities. 


EXPERIMENTAL RESULTS 

Typical still photographs of burning magnesi-um ribbons are shown 
in figure 3. Figure 3(a) shows the ribbon twisted horizontally through 
the flame. Figure 3(b) is a medium exposure showing the ribbon less 
clearly than figure 3(a), but giving more en^hasis to the flame struc- 
ture. Figure 3(c) is a much heavier exposure; the igniter and the un- 
burned ribbon are visible. Magnesium oxide smoke is seen streaming up 
from the fire, but the structural details of the flame are obliterated. 
The intense inner regions surrounding the ribbon in both figures 3(a) 
and (b) are surrounded by regions of somewhat lower luminosity which in- 
clude, although only under the ribbon, a very thin region of high in- 
tensity. Figure 4 contains a schematic drawing indicating the various 
regions and shows microphotometer traces made across the paths indicated. 
The microphotometer traces indicate that, while the outer zone a is 
not an optical illusion, its intensity is not as high as that of the 
inner zone c. 

The oxide ash left after a fire usually supported itself and spanned 
the gap. It maintained the general structure of the original ribbon, 
although it was fluffier for shorter burning times. For combustion in 
ordinary air, about 18 percent of the original magnesium remained in the 
ash as the oxide. The ash from the igniter thread also tended to pre- 
serve the structure of the original metal. 

The experimental burning times are presented in figure 5. The time 
decreases with increasing oxygen concentration for all three diluents. 

The decrease is somewhat less for the nitrogen-inerted case than for the 
argon case, and is appreciably less for helium than for either argon or 
nitrogen. An indication of the precision of the data is given by verti- 
cal lines indicating the spread of the two to four individual tests 
which are included in each average value plotted as a datum point. 
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The effect of moisture in the s\jrrounding atmosphere on h\nrning 
time is indicated in figure 6. Amounts of water vapor, up to saturation 
at room tenperature, were added to the mixtures of 21.5 and 43.7 percent 
oxygen in argon. The burning time decreases with increasing amounts of 
water, although the decrease in burning time per unit of added moisture 
is greatest for small quantities of water vapor. The precision of the 
measurements for low-moisture contents is generally poorer than the pre- 
cision indicated in figure 5 for equivalent burning times. As a result, 
the faired c\jrves are somewhat less significant than might be hoped; the 
qualitative effects are clear. 


DISCUSSION 

One of the objectives of this investigation was to learn whether a 
vapor -phase or a surface reaction controls the process under examination. 
The following evidence seems to indicate that a gaseous reaction occurs: 

(1) Distinct zones are visible in the photographs, somewhat like 
those of liquid hydrocarbon diffusion flames. (The outer luminous zone 
a, fig. 4, has tentatively been attributed to the recalescence of mag- 
nesium oxide.) 

(2) The light intensity, as indicated on the microphotometer traces, 
was greater above and below the ribbon than right on it (fig. 4). 

(3) The ribbon did not seem actually to liquify, although it did 
soften and twist; it supported itself on medium spans (for exanple, 
across the approximately 2.3 cm between the igniter and the bracket), 
althovigh it sagged on somewhat broader spans. The ash retained the 
general shape of the original ribbon. While these properties might in 
part belong to the oxide deposited on the sirrface rather than to the 
magnesium alone, the observations sviggest that the temperature at the 
metal surface is in the range between the melting (923° K) and the boil- 
ing (1393° K) points of magnesium, well below observed flame temperatures 
(>3000° K). 

(4) An appreciable loss of magnesium from the ribbon occurred d\jring 
the burning process, the oxide ash showed no indication of fusion, and 
the ash was fluffier for shorter burning times (faster rates). These ob- 
servations suggest the evaporation of magnesium at a temperat-ure well 
below that of the melting point of the oxide (>2500° K) . 

(5) Photographs of igniter threads burning in oxygen without ribbons 
showed flames of the same general size and shape as did threads biorning 
under similar conditions in chlorine. To produce flames in chlorine, 
heat had to be s-upplied continuously to the igniter and its thread; the 
magnesium apparently melted. This observation suggests a metal volatili- 
zation process for the chlorine flame, and, hence, from the similarities 
of the photographs for the two cases, a vapor -phase mechanism for the 
formation of the oxide. 
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In view of the evidence for the occurrence of a vapor -phase proc- 
ess in the burning of magnesium^ an attempt was made to compute burning 
rates, and, hence, burning times, by a method that has proved highly 
successful for single liquid fuel droplets. The two-region "stagnant 
film" approach of Spalding and later of Graves was adopted (ref. 19 and 
2l). Preliminary calculations (neglecting dissociation) suggested ex- 
cessive temperatures (T = 5X10^ multiplied by the mole fraction of 
oxygen); however, experimental measurements suggested a flame tempera- 
ture of about 3100^ K (refs. 15 and 16), approximately the boiling point 
of magnesium oxide. Since the oxide is believed to dissociate almost 
completely upon vaporization (ref. 24), dissociation could not be treated 
by the technique of reference 21. Rather it was decided to abandon the 
concept of an infinitesimally thin flame front and to consider a finite 
reaction zone of high-temperature diffusion. The model is shown in cross 
section in figure 7. It consists of three concentric cylindrical zones, 
with boundaries designated A, B, B*, and C. The metal is vaporized from 
the ribbon surface A at about the boiling point of magnesium. Then 
three zones AB, BB^, and B'C are established: 

(1) A zone AB of radially increasing temperature; the metal vapor 

is heated and heat is conducted to the metal surface to 
vaporize more metal. 

( 2 ) A zone BB^, a constant temperature reaction zone at the boiling 

point of magnesium oxide. Interdiffusion of the oxygen and 
magnesium occurs. The result is a chemical equilibrium be- 
tween condensed oxide, magnesium vapor, and oxygen all at the 
boiling point of the magnesium oxide in a stagnant film of 
inert gas. 

(3) A zone B*C of radially decreasing temperature, consisting of a 

stagnant film of inert gas with oxygen diffusing inward and 
being heated and with condensed oxide moving out and being 
cooled. [In the treatment of the spherical case, a steady- 
state solution is obtained as the radius of the C boundary 
approaches infinity; however, for the treatment of the cylin- 
drical or flat -plate cases’, a finite boundary condition is 
required. ] 

The general approach and the significant assuiip)tions involved in 
the calculation will be indicated here. The symbols used and the details 
of the calculation are presented in appendixes A and B, respectively. 

The physical assuii 5 )tions are: 

( 1 ) The metal ribbon can be adequately treated with cylindrical 
geometry; end effects are neglected. 
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(2) The temperature does not exceed 3200^ K; this follows from the 
assun5)tion of complete dissociation of magnesium oxide upon vaporization. 

(3) The regions of highest liaminosity correspond to the reaction 
zone; specifically, the apparent flame diameter has heen taken as the 

boundary. Some indication of the effect of assuming other boundary 
conditions will be given. 

(4) The effects of free convection are neglected in setting up fuel- 
flow, oxygen -diffusion, and heat-transfer equations for the steady-state 
process. 

(5) The required thermal conductivities and diffusion coefficients 
may be calculated, as no experimental data exist for the high-temperature 
values required. The thermal conductivity of the gases in the AB zone 

is taken as that of magnesium vapor (Mg) and in the B’C zone as that of 
the appropriate inert gas. Diffusion in the BB' and B*C zones is con- 
sidered as the diffusion of oxygen into a stagnant film of the appropri- 
ate inert gas . 

(6) In the BT zone, as in reference 25, the film thicknesses for 
heat and mass transfer are assumed to be equal. 

(7) The pressure throughout the system is at 1 atmosphere. 

The calculation of burning time is carried out by computing the 
burning rate in terms of fuel consumed in moles per second; this calcu- 
lation requires simultaneous solution of a number of differential equa- 
tions for fuel-flow rate, oxygen diffusion, and heat transfer. These 
are organized zone by zone and integrated using the boundary conditions 
presented in the following table: 


Boundary 

Radius , 
cm 

Temperature , 
°K 

Partial pressure 
of oxygen, 
atm 

A 

r^(metal surface) 

1393 (b.p. Mg) 


B 

r^( eliminated) 

3200 (approx, b.p. MgO) 

bd 

II 

0 

B' 

r^ T ( luminous surface ) 

3200 (approx, b.p. MgO) 

t ( calc • ) 

C 

rQ(not used) 

300 (ambient) 

p^( ambient) 


The following equations, derived in appendix B in order to calcu- 
late burning times, are indicated here in the integrated forms finally 
en5)loyed in computation. 
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AB Zone . - The heat-transfer equation for the metal v^apor zone AB 


IS 


W« r-n 

^ in ^ = B 

2rt2 — 


( 1 ) 


where B is a definite integral over the teicperature range and involves 
the thermal conductivity and thermodynamic properties of the metal vapor 
in practice, a graphical Integration is performed. 

BJC_^one. - The heat-transfer equation for the B’C zone is 


% T 

In 

2nl r^i 


= C 


( 2 ) 


where C is a definite integral over the temperature range and involves 
the properties of oxygen, the combustion products, and the inert gas; 
again, a graphical integration. 


Also in B*C zone the oxygen diffusion equation is 

(1 - Pb.)] (5) 

Where F, a fianction of In (l - Pg,), is a definite integral of the 
diffusion eqiiation over the range of the partial pressxire of oxygen. 

When Wj. is eliminated from equations (2) and (3), In (l - p^,) 

may he determined for each condition of inert and ambient oxygen 
concentration. 

BB' Zone . - The oxygen diffusion equation in the BB' zone is 

gj in ^ = B’ [in (1 - pg,)] (4) 

B 

where B' is a definite integral of the diffusion eqtiation over the 
range of partial pressure of oxygen. When values of In (l - Pg,) ob- 
tained by combining eqviations (2) and (3) are \ised, B' may be evaluated 
in equation (4). 

From eq\iations (l), (2), and (4), a series of expressions as de- 
rived in equation (B12), (appendix B), maybe obtained, including 
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Wj B+B* B+B* 

2nl T ’'c 

in — in - 

A A 

Then 


t = 4.85X10-5 g /g, - (6) 

or any of a number of similar expressions. The mass per unit length of 
the magnesium ribbon is 0.0074 gram per centimeter and the molecular 
weight of magnesium is 24.3. The definite integrals over temperature 
in the AB and the B'C zones, B and C, respectively, are relatively small 
as compared with the definite integral over partial press\jre of oxygen 
B ' in the BB ' zone for oxygen concentrations above 40 to 50 percent . 
Since B, B', and C can be evalmted by equations (B5) , (B7), and 
(Bll) of appendix B (values are given in fig. 8), the burning time will 
be determined by the values assigned to the radii of the boundaries of 
the zones. The radius r^ is associated with the ribbon; in particular, 

the radius of a cylinder of equal surface area (0.1025 cm) has been 
adopted; an alternative is the use of the radius of equivalent cross 
section (0.0386 cm), and the difference in the result is an increase in 
con 5 )uted burning time t by a factor of about 2. The radius r^i has 
been taken as one -half the apparent flame diameter; the data are shown 
in figxire 9. Val-ues have been taken as independent of oxygen concentra- 
tion; in any case, the error of such an assumption will be insignificant 
on an absolute basis and of second-order effect on a relative basis. 
Moreover, it should be noted that while free convection has been neg- 
lected in the mathematical derivation, the position of this experimental 
boundary may be strongly affected by convection. 

Figure 10 shows the effects on the con 5 )uted bTorning time and zonal 
structure of the ass\m 5 >tion that r^, rather than I’g i ^ is properly rep- 
resented by a constant diameter. Again, the effect in t on an absolute 
basis is small especially for oxygen concentrations above 40 to 50 per- 
cent; relative values of t axe decreased to some extent in the lower 
ranges of oxygen concentration. 

AlthoTigh both experimental and computed burning times are averages , 
the averages are not taken in precisely the same way. There is reason 
to believe that the burning rate which is equivalent to the fuel flow 
Wf may be a function of r, and indeed eq'uation (5) shows that it will 

be, xmless the ratios of the diameters of the cylindrical zones remain 
constant. However, in the experiment a ribbon was used, and therefore 


m 


M2it2t 


( 5 ) 
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an assumption of constant siirface area and, hence, constant effectiv'e 
is plausible. In the calculation, the implicit assumption is that 

the cylinder is ignited over its entire length Instantaneously and burns 
uniformly. In the experiment, the burning time includes time for the 
flame to spread along the ribbon and takes little account of the fact 
that some portions are coii 5 )letely burned out in less than the burning 
time; a rough estimate is that these effects increase the observed times 
oi by a factor of about 2 which may also have a second-order dependence on 

oxygen concentration. 

The experimental 100-percent-oxygen value may better be considered 
as B.n i^per limit than as a true value because of the large calculated 
effects of traces of inert gas and/or the possibility that heat-transfer 
processes take control for very short times. The behavior of the com- 
puted curve for very high oxygen concentrations is shown by the dashed 
portion between 90 and 100 percent in figure ll(a). 

Figures ll(a), (b), and (c) show the results of carrying out the 
calculation outlined previously for each of the inert gases, argon, 
nitrogen, and helium, respectively. Each figure gives the appropriate 
experimental curve from figure 5, the calculated curve as obtained in- 
corporating the assunqptions and boundary conditions specified, and a 
normalized computed curve for comparison purposes. Normalization in- 
volved the multiplication of the computed curve by the ratio, (experi- 
mental value, 50 percent 02)/(computed value, 50 percent Og). The values 

of the ratio used are 2.13 for argon, 2.04 for nitrogen, and 2.72 for 
heliiam; the constant for the heliiim case depends significantly \: 5 )on the 
concentration selected for the normalization procedure. As pre- 
sented, no physical significance is attached to the ratio; however, 
factors of the magnitude involved could be obtained by decreasing r^, 

by decreasing the high-ten 5 )erature diffusion coefficient » , or by 
considering the spreading time of the flame. 

The correct prediction of the effect of oxygen concentration is the 
matter of major significance here; that is, the prediction of correct 
relative rates. In addition, calculation of the proper order of magni- 
tude for the actual time is gratifying; although the absolute agreement 

within a factor of 2.5 for argon and nitrogen must be deemed fortuitous. 

The argon curves (fig. 11(a)) show the best fit between experiment 
and calculation, followed quite closely by the nitrogen case (fig. 11(b)). 
The somewhat less satisfactory agreement for nitrogen is most readily 
attributed to the fact that nitrogen is not truly inert in this process; 
magnesium nitride is known to form, and traces of ammonia were present 
after runs. Nitrogen acting as reactant would tend to reduce burning 
time at low oxygen concentrations as is observed. However, the close 
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agreement between the nitrogen and argon experimental curves shows that 
here those effects have been eliminated which in high-temperature slow- 
oxidation experiments cause magnesium to oxidize more rapidly in air 
than in pure oxygen (ref. 18). No significant difference was obtained 
in the calculated curves for nitrogen and argon. Although the nitrogen 
curve is very slightly disposed in the direction indicated by the exper- 
iments, the differences are small compared with the assumptions involved, 
particularly in the calculation of the high-temperature diffusion 
coefficients. 

The helium curve (fig. 11(c)) shows considerable departure from 
experiment at low oxygen concentrations (these mixtures are also of 
relatively low density compared with the others considered). A con- 
siderable part of this departure is no doubt due to the difficulties of 
con 5 )uting thermal conductivities and diffusion constants for mixtures 
of constituents varying widely in molecular weight. Efforts to account 
for concentration effects on the properties involved are of doubtful 
validity and, in any case, show little improvement in the final results. 

The calculations with the inert gases do predict that helium should 
show shorter burning times than argon and nitrogen and this prediction 
is in agreement with the experiments. In comparison, the effect of in- 
ert gases in reducing the burning velocities of homogeneous hydrocarbon 
flames appears as Ng > A > He (ref. 23); from this investigation, if 

there is an analogy between burning velocity and the reciprocal of the 
burning time, Ng - A > He. However, here it is known that nitrogen is 

not con 5 )letely inert and would tend to produce a shift in the direction 
observed; from Ng > A for the hydrocarbon case to Ng = A for the 

magnesium case. Furthermore in the homogeneous case, flame ten 5 )eratures 
are not equal as has been assumed herein. 

The effect of water is less satisfactorily explained in terms of 
the thermal and diff-usional model presented. It would appear more prob- 
able that moisture affects the chemical or diffusion processes at either 
B or B* or at both boundaries than that it changes the gross proper- 
ties within an entire zone. Water vapor produces an increase in the 
oxidation rate in these experiments as it did in the high-temperature 
slow-oxidation-rate studies on magnesium (ref. 18). 

Nevertheless, in spite of certain anomalies between the experi- 
mental data and the coinputed burning times, the general agreement 
strongly supports the basic assunp)tion of a vapor -phase reaction. 


SUMMARY OF RESULTS 

Magnesium ribbons were burned in mixt\ires of oxygen, 17 to 100 per- 
cent by volme, with argon, nitrogen, helium, and argon - water vapor. 
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Experimental burning times were determined, and some details of the 
physical mechanism of the combustion process have been investigated. 

The following results were obtained: 

1. The burning time decreased sharply with increasing oxygen 
concentration. 

2. The burning time in argon-oxygen mixtures decreased in the pres- 
ence of moisture, with traces of water vapor showing the greatest effect. 

3. A vapor-phase mechanism for the combustion of magnesium ribbons 
was very strongly indicated. 

4. Stagnant -film concepts involving a reaction zone of finite 
thickness permitted the calculation of relative burning time in terms 
of the properties of the oxidizing atmosphere. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 27, 1954 
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APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 
A, B, B', C boundaries of zones, fig. 7 
AB, BB', B'C zones between appropriate boundaries 


B 


B' 


D 

F 


H 


definite integral over temperature in AB zone, eq. (B5), 
moles/ ( cm) ( sec ) 

definite integral over partial pressure of oxygen in BB* 
zone, eg. (Bll) , moles/(cm) (sec) 

definite integral over temperature in B*C zone, eg. (B7), 
moles/ (cm) (sec) 

diffusion coefficient of oxygen, sg cm/sec 

definite integral over partial pressure of oxygen in B’C 
zone, eg. (BIO) , moles/(cm) (sec) 

fuel flow per unit area, moles/(sg cm) (sec) 

heat required for process per mole of fuel, cal 

enthalpy of fuel, product, oxygen at T, cal/mole 


^,A 


\,C> ^02,C 

k 

L 

I 

M 

m 

n 

P 


enthalpy of fuel at = 1393^ K, cal/mole 
enthalpy of product, oxygen at Tq = 300^ K, cal/mole 

thermal conductivity, cal/(sg cm) ( sec) (^C/ cm) 

heat of vaporization of fuel at Tq = 1393° K, cal/mole 

length of cylinder, cm 

molecular weight of magnesium, 24.3 g/mole 
mass of magnesium rihhon, g 

number of moles per mole of fuel; stoichiometric relation 
total pressure, atm 
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p 

partial pressure of oxygen^ atm 

Q 

heat of combustion of fuel at Tq = 300^ K, cal/mole 

R 

universal gas constant, 82.05 (cu cm) ( atm)/ (^C) (mole) 

r 

radius , cm 

£ ^ 

absolute temperature, 

t 

con 5 )uted burning time, sec 


fuel flow, moles/ sec 

Wn 

^2 

oxygen flow, moles/sec 
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APPENDIX B 


PREDICTION OF BURNING TIMES FOR MAGNESIUM RIBBONS 

The stagnant -film approach of Spalding (ref. 19) and Graves (ref. 

21) is modified to include a finite reaction zone of high -temperature 
diffusion in place of an infinitesimally thin flame front. Three sym- 
metrical, concentric zones surround a metal cylinder (fig. 7), the re- 
sult of the assumption of no convection. The metal is ass'omed to reach 
the boiling point during ignition. The temperature of the reaction zone 
is assumed not to exceed 3200^ K, approximately the boiling point of 
magnesium oxide; this assumption is based on the fact that magnesium 
oxide is believed to dissociate upon vaporization (ref. 24), and means 
essentially that the heat of reaction at the boiling point is equal to 
the negative of the heat of vaporization. The concentration of oxygen 
is assumed to be zero at the inner boundary of the reaction zone. The 
diffusion of magnesium vapor is not considered explicitly; it is assumed 
that the diffusion of oxygen is into a stagnant film of gas having the 
properties of the appropriate inert gas with a total pressure of 1 atmos- 
phere throughout the system. 

The burning time is calculated in terms of the fuel flow (burning 
rate) W^; the differential equations are set up zone by zone and inte- 
grated using the boundary conditions indicated in the text. The thermo- 
dynamic data used were as follows: heat of combustion and heat of vapori- 

zation of magnesium (ref. 26), enthalpy of oxygen (ref. 27), enthalpy of 
solid magnesium oxide from heat-capacity data (ref. 28), and enthalpy of 

5 

magnesium vapor from a heat capacity of — R. Diffusion constants were 

calculated by the methods of Hirschf elder , Curtiss, and Bird (ref. 29) 
and are presented in the following table: 


Inert 

Db3,(3200° K), 

Dg,c(300° K), 

k-n t pT 1 

(300° K), 

gas 

sq cm/ sec 

sq cm/ sec 

cal/cu cm 

Argon 

10.5 

0.196 

6.4x10“2 

Nitrogen 

10.9 

.207 

6.8 

Helium 

38.2 

.75 

15.2 


The thermal conductivities of the inert gases were also computed accord- 
ing to reference 29. For the thermal conductivity of magnesium vapor, 
values of the viscosity were obtained by the use of the Sutherland formu- 
la, with a Sutherland constant (2130) taken equal to 1.54 times the boil- 
ing point of magnesium and a diameter of 3.2X10-Q centimeters; thermal 
conductivities were obtained from viscosities by the Eucken relation 
(ref. 30). Figure 12 shows the values of the thermal conductivities 
employed. 
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Oi 

CO 

-4 


Consider the rates at which fuel and oxygen 


steady state hy the 


reaction Mg 


+ I Oj 


MgO: 


are consumed in the 


Gf 


Wf _ -^02 

2-nrl 2«r7nn 
^2 


(Bl) 


The negative sign signifies that fuel and oxygen are transferred in oppo- 
site directions. 


The equation for the diffusion of a single gas through a stagnant 
film (ref. 3l) 


to 

I 

o 

o 


^02 DP 1 dp 

2nrl “ RT (P - p) dr 


may be applied In the various zones as follows : 


Zone AB; 


^ = 0 and 
dr 



0 


(B2) 


Zone BB ' : 

D = Bggt the diffusion constant for oxygen into the appropri- 
ate inert gas 

Zone B’C: 


D = Db,c diffusion constant for oxygen into the appropri- 

ate inert gas 

The heat -transfer equation 

= - 2«rlk ^ (B3) 


where H is the heat requirement per mole of fuel, may be applied in 
various zones as follows: 
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Zone AB: 

H = li 4- ii£» — 

k = k^ the thermal conduct ivity of magnesium vapor, Mg 
is negative since heat is conducted inwardly. 

Zone BB ’ : 

dT/dr is zero since the zone is assumed to he at the boiling 
point of magnesium oxide; H is also equal to zero. 

Zone B^C: 

H = Q - Dp (hp - hp^Q) + 

k = kg.c the thermal conductivity of the appropriate inert 
gas (A, Ng, He) 

HWf is positive since heat is conducted outwardly. 

The details of the application of equations (Bl) , (B2) , and (B3) 
are presented for each zone as follows: 

AB Zone . - The heat -transfer equation in zone AB is 

[l + hf - hf , a] Wf = ^ab ^ (B4) 

The term in brackets is the heat which must be conducted per mole in- 
ward from the reaction zone to vaporize the metal and raise the metal 
vapor to the temperature T, assuming the metal reaches its boiling point 
during ignition. The area of the cylinder is 2jtrl. The heat effects 
involved in the passage of some magnesium oxide to the ribbon in the ex- 
periment have been neglected. 

Equation (B4) upon rearranging, applying the boundary conditions 
at A and B, and integrating yields 


pr- 


B 


^f T ^B 

In 

2nl 


Wf 

2rtrl 


dr = 


■Tg = 3200 


*^AB 


= 1393 


[L + hf. - hf^A] 


dT H B 


(B5) 
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The graphical integration used here to calculate B is independent of 
the properties of the s\arrounding atmosphere. 

B ' C Zone . - The heat -transfer equation in zone B'C is 

[q - n^ih^ - hp^c) + - ^02,c)]Wf = - kB-c 2«rZ § (B6) 

The term in brackets is the heat which must he transferred outward from 
the reaction zone in the combustion of 1 mole of fuel. The dissociation 
of oxygen (about 10 percent at 3200° K) is neglected. All magnesium 
oxide is assumed to pass outward through B’C zone for this heat balance. 

Equation (B6) upon rearranging, applying the boundary conditions at 
B' and C, and integrating giyes 

^Tg, = 3200 

^5:^ 4TS C 

T = 300 + »02(»02 - l>02,c)] 

(B7) 

The definite integral C may be obtained graphically but has not been 
used explicitly because of physical assumption (3) in the text which 
identifies the apparent flame diameter with the'^' boundary. 

The equation for the diffiision of oxygen through the stagnant film 
in the B’C zone is 


2nl r^, 


^f ^ _ ^^2 ^ ^B’C^ 1 ^ 

2jtrZ “ " 2«rZnQ “ tiq RT (P - p) dr 
2 2 


(B8) 


Elimination of from equations (B6) and (B8)^ which implicitly 

assumes that the film thicknesses for heat and mass tranfer are equal, 
produces 


■^’C ^ 0 ^ 

[q - "p(kp - + ’^Og^^Og ■ ^Og^c^] ^ ^ ^ 

Since k and d/t have approximately the same tenperature dependence, 
M/D may be considered a constant (table in this section). Equation 
(B9) upon integrating, applying the appropriate boundary conditions at 
B’ and C, and inposing the condition that P = 1 atmosphere yields 
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PTg, = 3200 


dT 


- n (h _) + n- (h_ ^ 

P P,C' Oo' Oo 


Tq = 300 


^0p,C^J 


^B'C^ 1 ^ ^ ■ PB' _ ^ 


^B’C ’^02^'^ ^ 


(BIO) 


Evalmting the left side graphically gives values of In (l - pg,) in 
terms of the properties of the surrounding atmosphere. 


BB' Zone. - The eqmtion for oxygen diffusion in the BB' zone upon 
integration and application of the boundary conditions at B and B' 
(including Pg = O) and the condition P = 1 yields 


W r 

f , ^ 

2 «l t-q 


nn RT 


In (1 - pg , ) = B ’ 


(Bll) 


The definite integral B' may be evalmted directly with values of 
In (l - Pb>) obtained from equation (BIO). The temperature throughout 

the BB’ zone is taken as 3200° K. Values for are taken as those 

for oxygen diffusing into the appropriate inert gas; the assvmiption is 
that the effects of magnesium vapor will be negligible or, at least, 
about the same in every case. 

Then, from equations (B5) , (B7), and (Bll), a series of expressions 
for Wj may be obtained; 




2rtl 


B B' 




C B+B' B + B’+C 



B’ + C 



(B12) 


In addition to the relatively simple treatment presented, the more 
con^jlex equations (ref. 3l) for the diffusion of two vapors (here mag- 
nesiimi and oxygen through the inert gas) through a stagnant film were 
tentatively considered for the reaction zone. In this procedure, a par- 
tial pressure for magnesium of 0.9 atmosphere at the B boundary was 
assumed; this implies a temperature at the metal surface slightly below 
the boiling point of magnesium. The results of preliminary calculations 
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indicate decreased values for B', although the dependence of on 

oxygen concentration remains about the same. The over-all effect appears 
to he a somewhat decreased oxygen dependence in the regions of oxygen con- 
centration below about 50 to 60 percent^ a somewhat increased dependence 
upon oxygen concentration in the regions above about 80 percent , and a 
general increase in the calculated burning times. Such a solution still 
involves certain approximations and increases the complexity of the cal- 
culations appreciably without assurance of increased accuracy. A com- 
plete solution of the problem, including magnesium-vapor diffusion within 
the AB zone, is far more complex and is not warranted by the available 
data. 
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Figure 1. - Magnesium ribbon burning apparatus. 
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Figure 2. - Selected frames from motion pictures of burning magnesium ribbon. 
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(a) Short. 



(c) Long. 

Figure 3. - Burning magnesium ribbon, increasing exposures top to bottom. 
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Figure 4. - Schematic drawing of flame; microphotometer traces. 
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Figure 6. - Variation of turning time vith water concentration. 
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Flgiire 7. - Cross section of cylindrical stagnant-film, model with finite 
reaction zone; p, partial pressure. 


3277 


3277 



-4 


o 

0 

0 





W1 


Figure 8. - Yalues of definite Integrals defined in equations (B5), (B7), and (Bll). B, 0.48x10"^ mole per 

centimeter per second; C for argon, nitrogen, and helium: 0.30, 0.33, and 2.52x10“^ moles per centimeter 
per second, respectively. 
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Figure 9. - Variation of apparent flame diameter with oxygen 
concentration and inert gas. 


Radii of zones, cm Burning time, sec 


MCA TN 3332 


33 






Oxygen, percent by volume in surrounding atmosphere 
(a) Tjg, held constant. (b) r^ held constant. 

Figure 10. - Effect of boundary conditions on zonal structure and burning time. 



(a) Argon. 

Figure 11. - Variation of burning time with oxygen concentration and Inert gases. 
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Burning time, sec 



(b) Nitrogen. 

Plgure 11. - Continued. Variation of burning time with oxygen concentration and Inert gases. 
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(c) Helium. 

Figure 11. - Concluded. Variation of burning time with oxygen concentration and Inert gases. 
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